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ABSTRACT 


An investigation of the effect of basis selection on geometric nonlinear response prediction using 
a reduced-order nonlinear modal simulation is presented. The accuracy is dictated by the 
selection of the basis used to determine the nonlinear modal stiffness. The scope of this 
investigation is limited to structures which do not exhibit linear bending-membrane coupling, but 
do exhibit nonlinear bending-membrane coupling. This study considers a suite of available bases 
including bending modes only, bending and membrane modes, coupled bending and companion 
modes, and uncoupled bending and companion modes. Companion modes represent an 
alternative to membrane modes and capture some of the membrane behavior resulting from 
bending-membrane coupling. The nonlinear modal simulation presented is broadly applicable 
and is demonstrated for nonlinear quasi-static and random acoustic response of flat beam and 
plate structures with isotropic material properties. Reduced-order analysis predictions are 
compared with those made using a numerical simulation in physical degrees-of-freedom to 
quantify the error associated with the selected modal bases. Bending and membrane responses 
are separately presented to help differentiate the bases. 

1. INTRODUCTION 

The design of advanced aerospace vehicle components capable of withstanding high 
vibroacoustic environments is hampered by a lack of accurate and computationally fast methods. 
Such methods are required in the design phase to quickly assess the impact of design changes on 
high-cycle fatigue life. Linear analysis methods are often inappropriate as structures may 
respond in a geometrically nonlinear fashion. Therefore, the use of a nonlinear analysis is 
required. Complicated structural geometries dictate the use of a finite element analysis (FEA). 
The traditional FEA employing numerical simulation in physical degrees-of-freedom (DoFs), 
however, is computationally intensive and considered impractical in design environments where 
rapid prototyping is needed. For stochastic response, the computational burden is exacerbated by 
the need to perform probabilistic analysis, such as Monte Carlo simulation [1,2], which requires 
that multiple realizations of the response be computed to generate meaningful statistics. 
Therefore, alternative methods are sought which retain the level of accuracy required, yet are 
computationally efficient. 
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Reduced-order nonlinear finite element analysis methods gain their computational advantage 
over direct numerical simulation by transfonning the equations of motion in physical DoFs to 
modal coordinates. Consequently, system size is significantly reduced and can be solved in a 
time efficient manner by means of numerical simulation or by an equivalent linearization 
approach, as dictated by the desired fidelity of the analysis. In the recent past, equivalent 
linearization procedures [3, 4] have been shown to be applicable [5, 6] to this class of problems, 
albeit in an approximate sense. This study focuses on reduced-order numerical integration 
analysis. 

The problems of interest in this paper are those in which the structure responds to imposed loads 
in a geometrically nonlinear (large deflection) static and random fashion. This nonlinearity is 
due to bending-membrane coupling and gives rise to membrane stretching when out-of-plane 
loading is applied. Consideration of structures which exhibit linear bending-membrane coupling 
is outside the scope of this work. Structures which exhibit linear bending-membrane coupling 
include curved structures and non-symmetrically laminated composites with a non-zero coupling 
stiffness matrix [f?] . In the recent past, a significant amount of research has been perfonned in 

reduced-order methods development that is applicable to the class of problems of interest. These 
methods may be viewed as being in one of two categories; those in which the nonlinear modal 
stiffness is directly evaluated from the nonlinear finite element stiffness matrix (so-called direct 
methods), and those in which the nonlinear modal stiffness is indirectly evaluated. Direct 
methods are typically implemented in special purpose finite element codes in which the 
nonlinear stiffness is known. The work by Mei et al [7] for flat configurations and the expansion 
of that approach for curved panels by Przekop et al [8, 9] are good examples of the direct 
method. In physical coordinates, the membrane displacements may be statically condensed into 
the bending displacements by neglecting in-plane inertia, and thus eliminate the need for 
membrane modes in the modal basis. This approach is sometimes referred to as “direct physical 
condensation” [10]. Alternatively, modal in-plane inertia may be neglected once the system of 
equations is transformed to modal coordinates, in an approach subsequently referred to as “direct 
modal condensation.” In this variation, the membrane behavior must still be represented in the 
basis. 
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For the analysis of complicated structures, commercial finite element codes are often required 
due to their support for a large number of element types. Unfortunately, the nonlinear stiffness is 
not typically available, making implementation of the direct method not possible. The only 
known implementation of a direct method into a commercial code is due to Bathe and Gracewski 
[11], in the ADINA finite element program. Indirect stiffness evaluation methods arose from the 
desire to implement reduced-order nonlinear analyses within the context of any commercial 
finite element analysis. Examples of indirect stiffness evaluation approaches may be found in 
the work of McEwan et al [12, 13] and Muravyov and Rizzi [14]. The approach taken by 
McEwan involves solution of a series of nonlinear static problems, typically through application 
of static forces with distributions corresponding to bending modes only. Following a modal 
coordinate transformation, the nonlinear modal stiffness is determined by curve fitting the modal 
force - displacement relation. Because the nonlinear static problem uses prescribed forces, the 
displacement response is influenced by the effect of membrane stiffness resulting from nonlinear 
bending-membrane coupling. The effect of membrane modes is thus implicitly condensed into 
the bending modes. This is advantageous in that explicit inclusion of membrane modes in the 
modal basis is eliminated. The implicit condensation of membrane modes however renders the 
approach unable to directly determine membrane displacements. Recently a post-processing 
procedure involving a mapping technique, called an “estimated expansion basis” [15], was 
developed to mitigate this problem. Nevertheless, the implicit condensation method is capable 
of producing only cubic nonlinear modal stiffness terms [15]. Its general applicability to non- 
planar structures, in which the effects of quadratic stiffness and membrane inertia may be 
significant, has yet to be explored. 

By contrast, the approach undertaken by the authors solves a series of simple algebraic equations 
obtained from static nonlinear analyses using prescribed displacements obtained from a 
combination of basis vectors. For the problems under consideration, the low-frequency bending 
modes obtained from the linear eigenvalue problem are uncoupled from the high-frequency 
membrane modes. In past works by the authors [5, 14], the basis vectors were formed from only 
low-frequency bending modes. However, to more accurately represent the effect of nonlinear 
bending-membrane coupling, the membrane response must be explicitly represented through 
inclusion of some form of membrane response in the basis. Inclusion of membrane modes may 
be cumbersome because the task of identifying a particular high-frequency membrane mode 
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amongst a great number of computed modes is labor intensive. In recognition of that, Hollkamp 
et al [10, 15] and Mignolet et al [16] have developed “companion” or “dual” modes to help 
capture the effect of membrane response, without explicitly including membrane modes in the 
basis. 

For both direct and indirect stiffness evaluation approaches, the accuracy of the solution depends 
on the selection of the modal basis, through which the nonlinear modal stiffness may be 
determined. If an insufficient basis is selected, then the predicted dynamic response of the 
reduced-order model may significantly differ from that of the full nonlinear model. Thus, there 
is always a need to assess the appropriateness of the selected basis via comparison of the 
predicted reduced-order response with something other than a reduced-order method, e.g. 
experimental data or numerical simulation in physical DoFs. Recent work compared various 
reduced-order approaches with experimental data [10]. In this study, comparisons are made with 
numerical simulation in physical DoFs to permit identical specification of boundary conditions. 

This paper assesses the effect of basis selection on the response obtained from a nonlinear modal 
simulation, utilizing the authors’ indirect stiffness evaluation method. A suite of bases is 
considered including bending modes only, bending and membrane modes, coupled bending and 
companion modes, and uncoupled bending and companion modes. The effect of basis selection 
on the modal stiffness coefficients themselves is first investigated. Then, using these 
coefficients, the nonlinear quasi-static and random response of simple planar aluminum beam 
and plate structures under spatially uniform excitation is considered. These structures were 
selected to help keep the cost of the comparative physical DoFs simulation reasonable, yet retain 
the nonlinear bending-membrane coupling behavior of interest. The error associated with the 
modal basis selection is quantified for both the displacement and stress response. Bending and 
membrane responses are separately presented to help differentiate the bases. 

2. NONLINEAR MODAL SIMULATION 

The nonlinear modal simulation analysis employed in this work consists of several parts. One or 
more methods, to be discussed, are first used to obtain a modal basis. Following a 
transformation of the nonlinear system to modal coordinates, the modal stiffness coefficients are 
evaluated and the resulting coupled system of equations is numerically integrated to obtain the 
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modal displacement time history. These are transformed back to physical coordinates for post- 
processing, including stress recovery. 


2.1. MODAL COORDINATE TRANSFORMATION 

The equations of motion of the nonlinear system in physical DoFs may be written as 

MX(t) + CX(t) + F NL (X(t)) = F(t ) (1) 

where M and C are the mass and proportional damping matrices, respectively, X is the 
displacement response vector and F is the force excitation vector. As written, the nonlinear 
restoring force vector, F NL , contains the linear force KX and nonlinear forces, where K is the 
linear stiffness. 

A set of coupled modal equations with reduced DoFs is first obtained by applying the modal 
coordinate transformation X - Oq to Equation (1), where q is the vector of modal coordinates. 
The modal basis matrix ® is typically formed from the eigenvectors obtained from Equation (1) 
using only the linear stiffness. For flat isotropic structures, these may include any combination 
of bending and membrane modes. In lieu of membrane modes, the modal basis may include 
“companions” related to the membrane response, as discussed in the next section. Generally, a 
small set ( L ) of basis vectors are included giving 

Mq{t) + Cq(t) + F nl (q l ( t ), q 2 ( t ), . . . , q, (0) = F(t) (2) 


where, for mass-normalized eigenvectors, 


M = <D r M® = [l] 

C = 0 T CO=[2C r coJ 

Fnl = ® T F nl 

f = o t f 


(3) 


and co r are the undamped natural frequencies and are the viscous damping factors. 
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2.2. MODAL BASIS SELECTION 


For the problems of interest in this paper, both bending and membrane behavior should be 
included in the basis selection since the large deflection nonlinearity couples their response. The 
basis vectors may be detennined via several methods. Bases corresponding to the bending and 
membrane response may be detennined through solution of the linear eigenvalue problem. Other 
basis vectors corresponding to the membrane response induced by bending-membrane coupling 
may be determined via alternative approaches, as discussed in Sec. 2.2.2. Lastly, while not the 
subject of this paper, basis vectors may also be determined via experiment or a hybrid scheme. 

2.2.1 Linear Eigenvectors 

Recall that for the problems of interest, the linear eigenvectors, obtained from Equation (1) using 
only the linear stiffness, are uncoupled and are either associated with low-frequency bending 
modes or high-frequency membrane modes. The selection of which bending modes to include 
depends on a number of factors including the excitation bandwidth, the spatial loading 
distribution, and even geometric and material properties. The selection of which membrane 
modes to include is less apparent than bending modes, as these high-frequency modes typically 
reside above the excitation bandwidth. Nevertheless, a reasonable starting point is to select the 
lowest membrane modes that are consistent with the spatial loading distribution and other 
physical properties. For example, for a uniform flat structure under uniform transverse loading, 
the membrane displacement will be anti-symmetric, so anti-symmetric membrane modes would 
be selected. Inclusion of both bending and membrane eigenvectors in the modal basis, either 
independently or in pairs, is subsequently referred to as the bending and membrane mode (BM) 
basis. Inclusion of only the bending eigenvectors will be referred to as the bending mode only 
(B) basis. In this study, the mass-normalized eigenvectors were obtained using 
MSC.NASTRAN normal modes analysis (solution 103). 

2.2.2 Companion Basis Vectors 

An alternative approach to using membrane modes is the use of so-called companion [10] or dual 
[16] modes. These modes are meant to represent the membrane behavior resulting from bending 
due to bending-membrane coupling. Previous authors utilized quasi-static approaches to 
determine the companion mode via nonlinear static analyses. For each companion mode, 
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Hollkamp et al [10] prescribed a displacement field corresponding to the bending mode of 
interest. Mignolet et al [16] applied two loadings, each having the distribution of a particular 
bending mode shape but with different magnitudes, to obtain the companion for that bending 
mode. In either case, the resulting modes contained only the membrane related behavior, but not 
the bending. The companion modes obtained need to be mass-normalized prior to their use. 

A new method is now presented for computing the companion mode using a dynamic analysis. 
An initial stress-free, out-of-plane perturbation of the mesh is first introduced to couple the 
bending-membrane response. The magnitude of the perturbation is chosen to be very small such 
that a normal modes analysis yields virtually the same bending eigenvalues and eigenvectors as 
that of the flat structure. An engineering rationale for specifying the shape of the perturbation is 
not evident. It was found that the companion modes change with varying perturbation shapes. 
In this investigation, the shape was chosen to be that of the first bending mode. The bending and 
high frequency membrane components remained identical to that of the flat structure, regardless 
of the shape of the perturbation. The MSC.NASTRAN normal modes analysis was again used to 
compute the mass-normalized eigenvectors, which now contain both the bending and membrane 
behaviors, but at the natural frequencies of the original flat structure. Direct inclusion of these 
eigenvectors in the modal basis is subsequently referred to as the coupled bending and 
companion mode (CBC) basis. 

A more consistent usage of companion modes, with respect to references [10, 16], is to separate 
the DoFs associated with the bending and membrane behaviors. In practice, since the bending 
behavior is unchanged, the original low-frequency bending modes are retained. The bending 
DoFs are set to zero in the newly obtained eigenvector to obtain the uncoupled dynamic 
companion. Since the uncoupled dynamic companion is essentially obtained by partitioning the 
mass-normalized CBC mode, it is no longer itself mass-normalized. Therefore, an additional 
step of mass-normalization is necessary. Inclusion of both the original bending and uncoupled 
dynamic companion modes in the modal basis, independently or in pairs, is subsequently 
referred to as the uncoupled bending and companion mode (UBC) basis. 
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2.2.3 Comparison of Membrane and Companion Modes 

Consider a clamped-clamped aluminum beam measuring 18-in. x 1-in. x 0.09-in (/ x w x h ) with 
the following material properties: 

E = 10.6 x 10 6 psi, G - 4.0 x 10 6 psi, p = 2.588 x 10^ 1 ^~ 

The beam was modeled in MSC.NASTRAN using 144 CBEAM elements. At the clamped ends, 
all DoFs are constrained. For the unifonnly distributed loadings to follow, the first six 
symmetric bending modes (eigenvectors 1, 3, 7, 10, 14 and 19) at natural frequencies of 58, 312, 
770, 1431, 2293 and 3354 Hz, respectively, were selected. The first six anti- symmetric 
membrane modes (eigenvectors 46, 81, 1 15, 153, 221 and 231) are at natural frequencies of 1 1.2, 
22.5, 33.7, 44.9, 56.1 and 67.3 kHz, respectively. The first four of these are plotted in Figure 1 - 
Figure 4. Also shown are the static companion modes obtained by Mignolef s approach [16] and 
the dynamic companion modes obtained by the method outlined above. It is clear that both static 
and dynamic companion modes significantly differ in shape relative to the membrane mode. The 
effect of their inclusion in the modal basis on the stiffness coefficients is next considered. 

2.3. INDIRECT STIFFNESS EVAFUATION METHOD 

The indirect stiffness evaluation method previously developed [14, 16] was used in this study. 
To summarize, the nonlinear force vector in Equation (2) may be written in the fonn 

L L L L L L 

FfjL (*7l J *?2’ • ' •’ 7i) ~ Z d j^j + Z Z ">,7/, + Z Z Z hjuVjMi r = l,2,...,L (4) 

j = 1 }=\ k= ] j = 1 k= ] l=k 

reducing the problem of determining the nonlinear stiffness from one in which a large set of 
simultaneous nonlinear equations must be solved to one involving simple algebraic relations. 
The linear, quadratic, and cubic nonlinear modal stiffness coefficients are written as dj , a jk and 

b jkl , respectively. 
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Figure 1: First anti- symmetric membrane and companion modes. 



Figure 2: Second anti- symmetric membrane and companion modes. 
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Figure 3: Third anti-symmetric membrane and companion modes. 



Figure 4: Fourth anti-symmetric membrane and companion modes. 
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For the prescribed static displacement fields [16] 


X* ^1 X c2 -^c3 +4Wl (5) 

the nonlinear forces evaluated using the MSC.NASTRAN nonlinear static solution (solution 106) 
are given as 


F Wi , (D + \_d\ J Qi + \_ a u J q$i + J *7i*7i*7i 

Fm, 2 = ® F M (-H) = -[,/i],i + [aii]^i-[hiii]^i (6) 

F M3 = (+<K^i ) = + \d [ ] q x + [a r n ] q x q i “ [&m] 

The first two of Equation (6) are used to evaluate \a r n ] . The first and the third of Equation (6) 
are then used to determine [d[] and [b' n , ] . It should be noted that the specified modal 
displacements in Equation (5) are single scalar quantities since the analysis is static. These 
modal displacements are specified such that the magnitude of the prescribed physical 
displacement field X c is physically meaningful. It has been shown that the stiffness coefficients 
obtained via this approach are not sensitive to the particular value of modal displacement 
specified [14]. The last of Equation (5) are obtained by specifying a scaled modal displacement 
such that q-aq. It was found in the course of this work that the stiffness coefficients obtained 

were also insensitive to the particular value of a used. In this paper, values of q = 7 x 1CT 5 and 
a = 1.25 were used for all modes in the basis. 

This procedure for evaluating the linear stiffness coefficients [d[ ] is advantageous over the 

earlier implementation [14] using the MSC.NASTRAN linear static solution (solution 101), as 
the linear solution does not include in-plane displacements, and therefore is not capable of 
handling thermal loads. The remaining nonlinear coefficients are evaluated through an additional 
series of nonlinear static solutions, as described in [14], 

2.3.1 Comparison of Modal Stiffness Coefficients 

For the clamped-clamped beam structure previously considered, the effect of including a 
particular set of modes in the modal basis is shown in Table 1 - Table 5. Table 1 shows the 
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diagonal terms of the linear stiffness coefficients for a set consisting of the first six symmetric 
bending modes ( d\ - d \ ) plus up to six membrane or companion modes ( cl] - d\; ). For each 

basis, the first six linear stiffness coefficients are equal to the eigenvalues of the selected bending 
modes. For the BM basis, the second six linear stiffness coefficients correspond to the 
eigenvalues of the selected membrane modes. The CBC basis lacks the second set of six 
coefficients as only six modes are included in this basis, i.e. six coupled bending and companion 
modes. It is interesting to note that the stiffness coefficients for companion modes of both UBC 
bases differ from each other and from the membrane. Thus, there is no physical meaning to the 
linear stiffness coefficients d ] - d\\\ associated with the companion modes for the UBC bases. 

Further, although comparable in magnitude to the membrane modes, the stiffness coefficients for 
companion modes of both UBC bases do not increase monotonically. For the CBC basis and 
both UBC bases, the off-diagonal terms are non-zero (not shown), indicating coupling of the 
linear stiffness where none should exist. 

Table 1 : Diagonal terms of beam linear stiffness obtained using five different modal bases. 

Shaded cells indicate bending eigenvalues. 



B 

BM 

CBC 

UBC 

ubc [16J 

(static) 

d\ 

1.318E+05 

1.318E+05 

1.318E+05 

1.318E+05 

1.318E+05 

d\ 

3.845E+06 

3.845E+06 

3.845E+06 

3.845E+06 

3.845E+06 

d\ 

2.341E+07 

2.341E+07 

2.341E+07 

2.341E+07 

2.341E+07 

d\ 

8.080E+07 

8.080E+07 

8.080E+07 

8.080E+07 

8.080E+07 

d\ 

2.075E+08 

2.075E+08 

2.075E+08 

2.075E+08 

2.075E+08 

d\ 

4.442E+08 

4.442E+08 

4.442E+08 

4.442E+08 

4.442E+08 

d] 

N/A 

4.990E+09 

N/A 

7.090E+10 

1.487E+10 

d\ 

N/A 

1.995E+10 

N/A 

7.621E+10 

7.520E+10 

dg 

N/A 

4.485E+10 

N/A 

2.587E+09 

5.846E+10 

4° 

N/A 

7.965E+10 

N/A 

2.057E+10 

7.803E+10 

d\\ 

N/A 

1.243E+1 1 

N/A 

N/A 

N/A 

i i 2 
a \2 

N/A 

1.786E+1 1 

N/A 

N/A 

N/A 


Unlike direct physical and direct modal condensation, or the implicit condensation in McEwan’s 
indirect approach, the effect of quadratic stiffness in the present stiffness evaluation approach is 
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not incorporated in the cubic stiffness, and must therefore be explicitly represented as indicated 
in Equation (4). Diagonal terms for the quadratic stiffness coefficients are shown in Table 2. 
The diagonal terms are expected to be zero for this flat structure, and essentially are for the B 
and BM bases. However, the CBC basis indicates significant coefficients for the first six modes 
( a \ j - a bb ). Stiffness coefficients for the companion modes of both UBC bases ( a 77 - a\l 10 ) are 

also significant, although their effect on the response is likely small for out-of-plane loadings. 
Interestingly, while the presence of anti-symmetric companion modes introduces significant 
diagonal terms for the CBC and UBC bases, the presence of anti- symmetric membrane modes 
does not affect the diagonal terms for the BM basis. The off-diagonal terms corresponding to the 
first (1 st symmetric bending) and seventh (1 st membrane or companion) modal equations are 
shown in Table 3, and indicate a coupling between modes for both UBC bases and the BM basis. 

Table 2: Diagonal terms of beam quadratic stiffness obtained using five different modal bases. 



B 

BM 

CBC 

UBC 

ubc [16J 

(static) 

a\ x q x 

-6.001E-1 1 

-6.653E-1 1 

6.937E+01 

-3.673E-1 1 

-1.737E-10 

2 

^22^2 

4.044E-10 

1.028E-10 

-1.474E+03 

1.874E-12 

4.471E-10 

^ 33^3 

-1.400E-09 

1.061E-09 

2.701E+03 

-5.210E-10 

7.226E-09 

^44^4 

-2.658E-09 

-1.834E-09 

3.879E+03 

2.047E-09 

-1.192E-07 

^55^5 

7.535E-09 

-6.914E-09 

-5.006E+03 

7.684E-10 

5.356E-09 


-3.101E-08 

-2.980E-08 

-6.295E+03 

-4.212E-08 

-2.676E-08 

^77*?7 

N/A 

8.864E-06 

N/A 

-7.245E-01 

6.642E+00 

8 

^ S 8*?8 

N/A 

-4.808E-05 

N/A 

1.767E+00 

1.570E+01 

9 

^99^9 

N/A 

-4.976E-05 

N/A 

-5.094E-01 

-7.283E+02 

10 

a io io^io 

N/A 

-1.408E-04 

N/A 

-5.227E-01 

-4.782E+03 

a l\ llt/ll 

N/A 

-4.697E-05 

N/A 

N/A 

N/A 

12 

a \2 12^12 

N/A 

5.307E-04 

N/A 

N/A 

N/A 


Relationships between quadratic modal stiffness coefficients were previously found in [14] and 
are given by 


1 — 9 7 

dy-j j j 

9 1 — 7 

dy-j 


( 7 ) 
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Equation (7) should be satisfied for only the significant nonlinear coefficients. For the quadratic 
terms presented in Table 3, the BM basis satisfies the first of Equation (7). For both UBC bases, 
however, the significant terms do not satisfy the first of Equation (7). The second of Equation 
(7) is not satisfied for any basis because these terms are not significant. 

Table 3: Selected off-diagonal quadratic stiffness terms obtained using five modal bases. 



B 

BM 

CBC 

UBC 

ubc [16J 

(static) 

a\ x q x 

-6.001E-1 1 

-6.653E-1 1 

6.937E+01 

-3.673E-1 1 

-1.737E-10 


N/A 

-2.320E+05 

N/A 

4.625E+01 

-2.398E+01 

a 11 q 1 

N/A 

-1.1 18E-06 

N/A 

3.312E-15 

3.608E-15 

a h<h 

N/A 

-1.160E+05 

N/A 

7.119E+10 

-8.107E+10 


N/A 

-3.946E-01 

N/A 

5.389E+01 

3.757E-01 

1 

& *1*1 *1 

N/A 

8.864E-06 

N/A 

-7.245E-01 

6.642E+00 


Table 4 shows the diagonal terms of the cubic stiffness coefficients for the same modal bases as 
above. For each basis, all coefficients corresponding to the bending modes (1-6) are identical. 
Like the quadratic stiffness, coefficients corresponding to the companion modes of both UBC 
bases are large relative to those corresponding to the membrane modes of the BM basis. 
Selected off-diagonal stiffness coefficients are shown in Table 5. The most significant difference 
between the non-zero stiffness terms for the BM and UBC bases is that the b \ n and the b 777 

terms are significant for the UBC bases and not for the BM basis. Conversely, the b\ nl term is 
significant for the BM basis but not for either UBC basis. 

Relationships between cubic modal stiffness coefficients were previously found in [14] and are 
given by 

T i _ >7 
°177 — °\ 17 

3 b] 17 = b\ in . (8) 

3 All =tfl 7 

Like the quadratic terms, Equation (8) should be preserved for only the significant nonlinear 
cubic coefficients. For the cubic terms presented in Table 5, the BM basis approximately 
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satisfies the first of Equation ( 8 ). For both UBC bases, however, neither the first nor the third of 
Equation ( 8 ) are satisfied since terms on only one side of the equation are significant. 

Table 4: Diagonal terms of beam cubic stiffness obtained using five different modal bases. 



B 

BM 

CBC 

UBC 

UBC [lbJ 

(static) 

b\\\ ( l\ 

5.067E+02 

5.067E+02 

5.067E+02 

5.067E+02 

5.067E+02 

^222 *7 2 

3.626E+04 

3.626E+04 

3.626E+04 

3.626E+04 

3.626E+04 

^333^3 

2.507E+05 

2.507E+05 

2.507E+05 

2.507E+05 

2.507E+05 

^ 444^4 

9.124E+05 

9.124E+05 

9.124E+05 

9.124E+05 

9.124E+05 

bid 

2.410E+06 

2.410E+06 

2.410E+06 

2.410E+06 

2.410E+06 

bid 

5.248E+06 

5.248E+06 

5.248E+06 

5.248E+06 

5.248E+06 

b-jrjQ-j 

N/A 

-4.897E-05 

N/A 

1.872E+02 

-2.340E+02 

^ 888^8 

N/A 

3.534E-04 

N/A 

-3.522E+00 

-1.121E+03 

^999 *?9 

N/A 

-3.199E-05 

N/A 

4.009E+01 

1.121E+04 

b\Q 10 10^10 

N/A 

5.931E-07 

N/A 

-1.499E+01 

-9.658E+04 

^1111 ll*7ll 

N/A 

2.756E-03 

N/A 

N/A 

N/A 

b] 2 12 12^12 

N/A 

-7.361E-03 

N/A 

N/A 

N/A 


Table 5: Selected off-diagonal cubic stiffness terms obtained using five different modal bases. 



B 

BM 

CBC 

UBC 

ubc [16J 

(static) 

b hd 

5.067E+02 

5.067E+02 

5.067E+02 

5.067E+02 

5.067E+02 

^117#l#7 

N/A 

-2.177E-06 

N/A 

-3.244E-1 1 

2.418E-09 

^177^1^7 

N/A 

-1.174E+03 

N/A 

-2.829E-06 

9.352E-07 

b-jj-jq-j 

N/A 

-6.761E-07 

N/A 

2.310E-14 

9.921E-15 

bud 

N/A 

7.559E-02 

N/A 

5.522E+03 

4.472E+04 

bwi^xq-] 

N/A 

-9.437E+02 

N/A 

8.816E+03 

-3.782E+03 


N/A 

-3.496E-01 

N/A 

5.390E-01 

1.671E+00 

b-m^! 

N/A 

-4.897E-05 

N/A 

1.872E+02 

-2.340E+02 


Further discussion regarding the choice of basis on the response prediction is reserved until 
Section 3. 
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2.4. NUMERICAL INTEGRATION AND ELEMENT STRESS RECOVERY 

Having the nonlinear force vector in Equation (4) fully defined, the coupled modal nonlinear 
equations of motion in Equation (2) may be solved. An equivalent linearization approach was 
presented in [5] to solve the system in an approximate manner. A more accurate approach 
utilized in this paper numerically integrates Equation (2) using a fourth-order Runge-Kutta 
method [17]. The resulting modal displacement time histories are transformed back to physical 
coordinates using the inverse modal transformation. 

Element stresses are recovered by post-processing the nodal physical DoFs directly within the 
finite element program. Since the stress post-processor of the finite element program is used, the 
stress calculation is identical to what would have been performed via a standard finite element 
analysis in physical DoFs. For a particular output time step, the element physical DoFs, obtained 
via the nonlinear modal simulation method, are applied to each element node as prescribed 
displacement fields in the MSC.NASTRAN nonlinear static solution and the element stresses are 
calculated. By repeating this operation for each output time step, the stress time history is 
determined. 

3. NONLINEAR STATIC AND RANDOM ACOUSTIC RESPONSE OF A BEAM 

The effect of modal basis selection on the nonlinear response is next considered for the same 
beam structure previously examined. It should be noted that the nonlinear modal simulation 
approach considered herein is applicable to arbitrary spatial and temporal loading distributions. 
Since the problems of interest originate from uniform acoustic loadings, the cases selected are 
limited to structures exposed to spatially unifonn loadings. Therefore, while the calculation of 
nonlinear stiffness coefficients as described above is independent of the loading, conclusions 
drawn about the modal basis selection, specifically the choice of which modes to include for a 
particular basis type, are load-specific. The merits of one basis type (B, BM, CBC and UBC) 
versus another, however, can be assessed independent of the loading distribution. For the 
uniform spatial distribution considered, two temporal variations are considered; a quasi-static 
loading and a band-limited random acoustic loading. Results from the nonlinear modal 
simulation are compared with those obtained via numerical simulation in physical coordinates. 
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3.1. NUMERICAL SIMULATION IN PHYSICAL COORDINATES 


Numerical simulation analyses in physical coordinates serves as the basis for comparison of 
results from nonlinear modal simulation analyses. The finite element program ABAQUS 
(version 6.4) was used to generate nonlinear displacement and stress time histories. The double 
precision explicit integration scheme with an adaptive time integration step (referred to as 
“element by element” in ABAQUS) was utilized for all analyses. The ABAQUS model 
consisted of 144 B21 elements. The choice of ABAQUS explicit analysis over alternative 
analysis methods, e.g. MSC.NASTRAN nonlinear transient solution (solution 129), was made 
based on its superior ability to simulate very long response histories. This choice introduced an 
inconsistency of element formulations used between the nonlinear modal and physical 
simulations. However, previous results (not shown) from a nonlinear static analysis in physical 
DoLs using MSC.NASTRAN (solution 106) and ABAQUS were nearly identical, indicating no 
significant difference between the different element formulations used. 

3.2. QUASI-STATIC RESPONSE 

In the limit of including all linear eigenvectors (the full set of BM modes) in the basis, results 
found using nonlinear modal simulation should be identical to those obtained via numerical 
simulation in physical coordinates because there is no modal truncation. The first purpose of the 
quasi-static response prediction is to investigate the effect of modal truncation. The second 
purpose is to investigate the spatial characteristics of results produced using each modal basis. 
Because the quasi-static response is considered, dynamic effects due to the modal mass and 
damping are not significant. Therefore, only the modal stiffness coefficients dictate the accuracy 
of the solution for a particular basis. Further, as these coefficients are the same for the quasi- 
static and subsequent dynamic analyses, the quasi-static analysis provides a benchmark against 
which the dynamic results may be judged. 

While it is possible to set the modal accelerations and velocities to zero in Equation (2) and solve 
for the modal displacements via Newton-Raphson (as per [14]), the more direct approach utilized 
here applies a quasi-static loading and numerically integrates to obtain the modal displacements. 
Once transfonned, the physical displacements and stresses obtained from the reduced-order 
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analysis are compared with those obtained from a nonlinear simulation in physical DoFs with the 
same quasi-static loading. 

3.2.1 Load Generation and Stress Post-Processing 

A uniformly distributed quasi-static pressure loading was applied to the beam to obtain the quasi- 
static response. The loading time history ramped from zero magnitude at time zero to its 
maximum level at 1.6384s, following which 0.5s of constant loading was applied. The response 
was obtained at 2. 1384s, after the decay of any transient behavior. 

The bending and membrane stress components were separated by averaging the upper and lower 
surface stresses to obtain the membrane component, and then subtracting that from the surface 
stress to obtain the bending component. This separation allowed greater insight into the effect of 
the various modal bases on the stress response. 


3.2.2 Results 

The effect of modal truncation was first considered. Three BM modal bases were used; a 4- 
mode basis consisting of the first two symmetric bending modes and first two anti-symmetric 
membrane modes (2+2), an 8-mode basis consisting of the first four symmetric bending modes 
and first four anti-symmetric membrane modes (4+4), and a 12-mode basis consisting of the first 
six symmetric bending modes and first six anti-symmetric membrane modes (6+6). 
Displacement results are presented at the % span (4.5-in. from the clamped end), and element 
stresses near the 14 span (4.4375-in. from the clamped end). These locations will help to 
elucidate the benefits and liabilities of the various modal bases under severe conditions as the 
membrane stress component is more significant there relative to the bending stress component 
than at other locations along the length, e.g., the clamped end or mid-span. 


The transverse displacement response error is shown in Figure 5 for a range of loadings. Results 
shown in this and in the subsequent three figures are relative to the physical DoFs solution from 
the ABAQUS explicit analysis, i.e., 


% Error = 


Physical - Reduced Order' 
Physical 


xlOO . 


(9) 
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Figure 5: Transverse displacement error at beam % span for three BM bases. 


It is seen that all bases except for 2+2 agree very well with the physical DoFs solution and 
asymptote to a small error. Since the subsequent random analysis in physical DoFs could only 
be perfonned using the ABAQUS explicit solution, the same solution was used here for 
consistency. Note that if the MSC.NASTRAN nonlinear static solution was used instead of the 
ABAQUS explicit analysis, a different set of small errors would be obtained. For example, for 
the 6+6 basis, at the highest level of 0.4608 psi, the error would be -0.07% using 
MSC.NASTRAN instead of +0.3%. Since the error is very small, Figure 5 should not be 
misinterpreted to indicate that the 4+4 basis is superior to the 6+6 basis at the highest loading 
level. Rather, the error in the 4+4 and 6+6 reduced-order analyses falls into an error band that is 
insignificant, regardless of the physical DoFs solution used for comparison. 


The membrane displacement response error is shown in Figure 6. Flere the 2+2 basis 
substantially under-predicts the response across the loading range, the 4+4 basis slightly over- 
predicts, and the 6+6 basis compares well across the range. The errors in bending and membrane 
stress components are shown in Figure 7 and Figure 8, respectively. In both cases, the 6+6 basis 
compares the most favorably across the loading range. The above results show that it is possible 
to achieve highly accurate displacement and stress response predictions using a truncated BM 
basis. 
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Figure 6: Membrane displacement error at beam % span for three BM bases. 
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Figure 7: Bending stress error at beam V* span for three BM bases. 


20 


10 I— 


□ 2 Bending + 2 Membrane 

A 4 Bending + 4 Membrane 
O 6 Bending + 6 Membrane 


LU 

t/> 

(/) 

o 


CO 

<D 

C 

0 

L- 

-Q 

E 

cu 


5 - 




□ 


o 

A 


□ 


o 


A 


□ 


o 


A 


-5 * ' I I 

10' 3 10' 2 10' 1 10° 

Pressure (psi) 

Figure 8: Membrane stress error at beam % span for three BM bases. 


Next considered are the effects of different modal basis variants having a comparable number of 
included modes. The following modal bases were considered; the first six symmetric bending 
and first six anti-symmetric membrane modes (BM 6+6), the first six symmetric bending modes 
(B 6), the first six coupled bending-companion modes (CBC 6), and the dynamic version of UBC 
having the first six symmetric bending modes and first four companion modes (6+4). A 
numerical instability for both UBC bases was encountered, limiting the number of companion 
modes for a particular excitation level. Radu et al [18] attributed this problem to the high 
frequency content of the companion modes. However, since the problem was not observed for 
the BM basis having comparable frequency content, a more likely reason has to do with the high 
degree of coupling between some bending and companion modes, as previously discussed. A 
static condensation in modal coordinates was proposed and demonstrated to mitigate this 
behavior [18], although no attempt was made to implement this scheme in the present work. The 
instability was severe enough, however, that results for the UBC static basis could not be 
obtained for the quasi-static loading case. For this reason, only the (6+4) UBC dynamic basis 
was considered. 


A single loading level (0.0576 psi) was applied and results are presented at locations across the 
beam semi-span to show how accuracy varies with location. Transverse displacements for each 
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modal basis are shown in Figure 9. The BM and UBC bases compare very well with the 
physical DoFs solution, while the B and CBC somewhat under-predict the response. By 
comparison, only the BM basis compares well with the membrane displacements shown in 
Figure 10. The bending-only basis B is incapable of predicting any membrane displacement, 
while the CBC basis correctly predicts the shape, but significantly under-predicts the magnitude. 
The UBC basis both over-predicts the magnitude and incorrectly predicts the shape. 



Figure 9: Quasi-static transverse displacements along beam semi-span for 4 modal bases. 

The bending and membrane stress behaviors are shown in Figure 1 1 and Figure 12, respectively. 
Bending stress results from all bases compare well with the physical DoFs solution. The nearly 
constant membrane stress is accurately captured by only the BM basis. The other bases predict 
incorrect magnitudes and spatial distributions. With regard to the bending-only basis, the 
membrane stress arises from the bending-membrane coupling in the stress recovery calculation, 
even though this basis predicts a zero membrane displacement. The fact that there are several 
points along the beam where membrane stress results from the B, CBC and UBC agree with the 
physical DoFs solution highlights the need to consider a number of spatial locations in assessing 
the accuracy of the modal basis. Note that only the BM basis accurately predicts bending and 
membrane displacement and stress response at all span-wise locations. 
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Figure 10: Quasi-static membrane displacements along beam semi-span for 4 modal bases. 



Figure 1 1 : Quasi-static bending stress along beam semi-span for 4 modal bases. 
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Figure 12: Quasi-static membrane stress along beam semi-span for 4 modal bases. 


A compact measure of an error along the span of the beam may be found from 


% Error = 


1 -Ai | ( Physical; - Reduced Order ) 
N y Physical f 


x 100 


( 10 ) 


where A is a number of points along the beam span used for the estimation. Figure 13 and 
Figure 14 present the transverse and membrane displacement errors, respectively. It is seen that 
BM basis provide the best reduced-order approximation through the entire range of excitation 
levels studied. From Figure 14 it can be also concluded that while the CBC basis outperforms 
the UBC basis, both incur unacceptably large errors. 

Figure 15 and Figure 16 present the bending and membrane stress errors, respectively. For the 
bending stress, the superiority of the BM basis is most evident at the highest loading. For the 
membrane stress, only the BM basis achieves errors less than 5%, while the B, CBC and UBC 
bases incur unacceptably large errors. 
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Figure 13: Quasi-static transverse displacement error integrated along the beam span. 



Figure 14: Quasi-static membrane displacement error integrated along the beam span. 
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Figure 15: Quasi-static bending stress error integrated along the beam span. 



Figure 16: Quasi-static membrane stress error integrated along the beam span. 

To summarize the findings from the quasi-static analysis, only the BM basis was able to capture 
the transverse and membrane displacements, and the bending and membrane stresses with high 
accuracy. The errors associated with the B, CBC, and UBC bases are generally smaller for the 
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transverse displacement and bending stress than they are for the membrane displacement and 
membrane stress. The errors associated with the membrane displacement and stress are so gross 
that the B, CBC and UBC bases are deemed unsuitable for their prediction for structures 
exhibiting nonlinear bending-membrane coupling. Consequently, if the membrane contribution 
is significant relative to that of bending, the B, CBC and UBC bases will also not be suitable. 
Therefore, some assessment of the membrane contribution must be made via an alternative 
analysis before blind application of the B, CBC or UBC bases. 

3.3. RANDOM ACOUSTIC RESPONSE 

For the nonlinear random response investigation, the beam was subjected to a unifonnly 
distributed acoustic loading with a bandwidth of 1500 Hz. Three loading levels were considered 
to span the response regime from essentially linear to highly nonlinear; a low overall sound 
pressure level of 128 dB (0.0072 psi RMS), a medium level of 146 dB (0.0576 psi RMS), and a 
high level of 164 dB (0.4608 psi RMS). Damping was chosen to be sufficiently high so that 
good comparisons could be made at the peaks of the PSD. A level of mass proportional damping 
was specified corresponding to 2.0% critical damping for the first symmetric bending mode. The 
damping specified was the same for both the reduced-order and physical DoFs simulation 
analyses. For the nonlinear modal simulation, a fixed time integration step of between 50 /is 
(for the lowest excitation levels) and 2 /is (for the highest excitation levels) was used. 

3.3.1 Foad Generation and Ensemble Averaging 

The same loading time history and ensemble averaging was used for the numerical simulation 
analysis in physical and modal coordinates. The loading time histories were generated by 
summing equal amplitude sine waves, each with random phase, within the specified bandwidth 
using a discrete inverse Fourier transform. This procedure was identical to that used in previous 
work [5] by the authors, so further details are omitted for brevity. The loading produced by this 
method has a Gaussian distribution. A sharp roll-off of the input spectrum practically eliminates 
excitation of the structure outside the frequency range of interest. 

For each load level, ten ensembles of displacement and stress response were generated lasting 
2.1384s each. The first 0.5s of each response record was discarded to remove the initial transient 
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response [5], resulting in response histories of 1.6384s in duration. For each simulation, the 
displacement and stress response was stored at every 50 jus , giving time records of 32,768 
points. A 32,768-point FFT was subsequently used to compute the power spectral density (PSD) 
function. 

3.3.2 Results 

The results shown are for the same bases considered in the quasi-static response analysis. The 
spectral results have been visually smoothed to more clearly observe the behavior. Transverse 
displacement PSDs are shown in Figure 17 for the medium loading of 146 dB. For this level, all 
modal bases agree well with the numerical simulation in physical DoFs, with the most significant 
differences at the third and fourth mode for the B, CBC and UBC bases. At this level, the 
transverse displacements obtained by the B and CBC bases are essentially the same. 



Figure 17: Quarter-span transverse displacement PSD of beam at 146 dB. 

The membrane displacement PSDs for this load level, shown in Figure 18, offer a different 
perspective. The BM basis is the only one to compare well with the physical DoFs simulation 
results. The CBC basis is inaccurate on two accounts; the magnitude is clearly incorrect and the 
shape of the frequency response mimics that of the bending response in Figure 17. Like the 
quasi-static response, the UBC basis over-predicts the response amplitude. Though the loading 
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spectrum falls off sharply at 1500 Hz, the membrane displacement response obtained using 
physical DoFs simulation and nonlinear modal simulations with B, BM, and UBC bases, extend 
to twice that frequency because of period doubling (not shown). The membrane displacement 
response computed using the CBC basis does not exhibit period doubling. 



Figure 18: Quarter-span membrane displacement PSD of beam at 146 dB. 

The bending stress PSDs at 146 dB are shown in Figure 19. The BM basis compares very well 
across the frequency range, even at the anti-resonances. The B and CBC stress PSDs are 
essentially the same, and under-predict the third and fourth resonances. The UBC basis is 
similar, but also under-predicts the fundamental. The membrane stress PSDs are shown in 
Figure 20. The BM basis compares remarkably well with the physical DoFs simulation. Again, 
the B and CBC bases predict nearly the same response to each other. Below about 630 Hz, the 
B, CBC, and UBC results fall on either side of the physical DoFs simulation, while above that 
frequency they all over-predict the response. 

At the highest loading of 164 dB, the transverse displacement PSD exhibits the peak spreading 
and shifting characteristic of spring hardening nonlinearity, see Figure 21. For this case, the 
bending-only and CBC bases indicate a higher degree of nonlinearity than the physical DoFs 
simulation results, as evidenced by the first and second peaks shifted to higher frequencies, 
greater broadening, and smaller magnitudes. Results from the UBC basis were not available due 
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to the aforementioned numerical integration problems. The BM basis results compare well 
across the frequency range, with some loss of peak response in the third peak. 



Figure 19: Quarter-span bending stress PSD of beam at 146 dB. 



Figure 20: Quarter-span membrane stress PSD of beam at 146 dB. 
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The membrane behavior shown in Figure 22 follows the observations made at 146 dB. The 
excellent agreement between BM basis and physical DoFs simulation results demonstrate the 
efficacy of the BM basis in the highly nonlinear dynamic response regime. 



Figure 2 1 : Quarter-span transverse displacement PSD of beam at 164 dB. 



Figure 22: Quarter-span membrane displacement PSD of beam at 164 dB. 

Bending and membrane stress PSDs at 164 dB are shown in Figure 23 and Figure 24, 
respectively. The BM basis most closely compares with the physical DoFs simulation bending 
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stress across the frequency range, with some loss in amplitude response in the peak near 850 Hz. 
The bending-only and CBC bases significantly under-predict the bending stress of the first and 
second modes and, above about 950 Hz, indicate different character than the BM basis and 
physical DoFs results. 



Figure 23: Quarter-span bending stress PSD of beam at 164 dB. 



Figure 24: Quarter-span membrane stress PSD of beam at 164 dB. 
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In Figure 24, the BM basis results compare favorably with the physical DoFs simulation results 
across the frequency range. The bending-only and CBC bases over-predict the membrane stress 
magnitude and width of the first peak. Above about 750 Hz, the character substantially differs 
from the physical DoFs results. 

To quantify the error in the root-mean-square of the response, the following measure was used: 


% Error = 


f RMS (Physical DoFs ) - RMS (Reduced Order) 
RMS (Physical DoFs ) 


xlOO . 


(ID 


The error in RMS response is an indicator of the error over the entire frequency range, and hence 
is not capable of representing the differences in the spectral character. Figure 25 and Figure 26 
present the error estimates for transverse and membrane displacements, respectively, and Figure 
27 and Figure 28 present the error estimates for bending and membrane stresses, respectively. 
While the RMS transverse displacement is estimated with least accuracy by the BM basis, it 
should be noted that this is primarily attributable to differences in the fundamental mode, which 
dominates the response. The 7% error in the RMS displacement transverse displacement is 
considered acceptable for a random analysis, particularly at the high nonlinear levels where the 
width of the confidence intervals is large [5]. 



Figure 25: Error in random acoustic RMS transverse displacement response at beam % span. 
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On the other hand, the RMS membrane displacement is captured most accurately by the BM 
basis. As shown in Figure 18 and Figure 22, the membrane displacement is dominated by the 
zero frequency component, which results from the fact that the membrane displacement 
oscillates between zero and some positive value. For the error in RMS membrane displacement, 
the CBC and UBC bases can be off by as much as one order of magnitude. 



Figure 26: Error in random acoustic RMS membrane displacement response at beam % span. 



Figure 27: Error in random acoustic RMS bending stress response at beam % span. 
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With regard to the stress response in Figure 27 and Figure 28, the BM basis provides the most 
accurate stress response, with a somewhat higher error in the bending stress at the medium load 
level. Again, this is likely a reflection of differences in the peak values. Bases B and CBC 
provide the lowest quality of stress estimation. For the bending stress, the RMS value is affected 
by several peaks of comparable magnitude, while the RMS membrane stress is dominated by the 
zero frequency component (see Figure 20 and Figure 24). 



Figure 28: Error in random acoustic RMS membrane stress response at beam l A span. 

4. NONLINEAR STATIC AND RANDOM ACOUSTIC RESPONSE OF A PLATE 

A plate was next considered to help understand the complexities of modal basis selection for a 
two-dimensional structure. Since the beam results presented in Sec. 3 indicated that the BM 
basis was superior over all the other basis variants considered, the nonlinear reduced-order 
results for the plate were computed using only the BM basis. 

The plate measured 14-in. x 10-in. x 0.04-in. (x x y x h ) with the aluminum material properties 
introduced in Sec. 2.2.3. The boundary conditions were simply supported on all sides. 
Specifically, along the 14-in. sides (in the x-direction), the transverse displacement DoF (w), the 
membrane displacement DoF in the y-direction (v), and the rotational DoF about the y-axis were 
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constrained. Along the 10-in. sides (in the v-direction), the w DoF, the membrane displacement 
DoF in the x-direction (u), and the rotational DoF about the x-axis were constrained. 

The ABAQUS model used for simulation in physical DoFs consisted of 8960 S4R elements 
measuring 0.125-in. x 0.125-in. The S4R element has one integration point, compared to four 
for the S4 element, and is required in the ABAQUS explicit analysis. The double precision 
explicit integration scheme with adaptive time step was used for both quasi-static and random 
acoustic response analyses. 

For the reduced-order analysis, the MSC.NASTRAN model used consisted of 8960 CQUAD4 
elements measuring 0.125-in. x 0.125-in. Up to the first twelve (eigenvectors 1, 4, 8, 11, 12, 19, 
22, 23, 28, 31, 35, and 40) doubly symmetric bending modes, i.e., symmetric in both x-y plan- 
form directions, were included in the basis. These occurred at natural frequencies 58, 216, 367, 
524, 531, 838, 985, 1003, 1140, 1308, 1452, and 1633 Hz, respectively. For the membrane 
modes, the selection was guided by the previous experience gained from the beam analysis, that 
is, membrane displacement components were either both anti- symmetric, or one was anti- 
symmetric and the other was zero. Up to the first twelve membrane modes were included 
(eigenvectors 426, 438, 618, 638, 778, 798, 894, 917, 935, 1142, 1152, and 1184). These 
occurred at natural frequencies of 15.2, 15.4, 21.2, 21.9, 26.1, 26.6, 29.6, 30.3, 30.8, 36.7, 36.9, 
and 37.9 kHz, respectively. 

4.1. QUASI-STATIC RESPONSE 

The effect of modal truncation on the quasi-static displacement results was first considered. The 
applied loading had a quasi-steady characteristic, as described in Sec. 3.2.1, and was uniformly 
distributed over the plate surface. Three BM modal bases were used; an 8-mode basis consisting 
of the first four symmetric bending modes and first four anti-symmetric membrane modes (4+4), 
a 16-mode basis consisting of the first eight symmetric bending modes and first eight anti- 
symmetric membrane modes (8+8), and a 24-mode basis consisting of the first twelve symmetric 
bending modes and first twelve anti-symmetric membrane modes (12+12). By analogy to the 
beam results, plate displacement results are presented at the % - % span location (3.5-in. from the 
boundary in the x-direction, and 2.5-in. from the boundary in the v-direction). 
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The transverse displacement (w) error and both membrane displacement (u in the x-direction, v 
in the v-direction) errors were computed as per Equation (9), and are shown in Table 6 for a 
range of loadings. It is seen that for the lowest quasi-static loading of 5.801 xlO’ 4 psi, the modal 
truncation has little effect on the displacement accuracy, and results computed with eight, 
sixteen, and twenty-four basis vectors are essentially the same. The center displacement at this 
level was w^/h = 0.016. As the quasi-static loading increases, the effect of nonlinearity 
increases, with a center displacement of w max //? = 0.408 for the 1.856 xlO'" psi load, and 
w max /fl =1-545 for the 0.297 psi load. The number of basis vectors used plays an increasingly 
important role as the load increases. The effect is more pronounced in the membrane 
displacement component errors than it is in the transverse displacement error. 


Table 6: Plate quasi-static displacement errors at (%, %) location, %. 



w 

u 

V 

5.801 xlO' 4 psi 

BM (12+12) 

-2.16 

-4.32 

-4.79 

BM (8+8) 

-2.17 

-4.39 

-0.18 

BM (4+4) 

-2.55 

-2.89 

-1.31 

1.856 xlO' 2 psi 

BM (12+12) 

-1.52 

-3.52 

-3.41 

BM (8+8) 

-1.41 

-2.55 

2.54 

BM (4+4) 

1.34 

4.23 

7.44 

0.297 psi 

BM (12+12) 

-0.48 

-3.89 

0.56 

BM (8+8) 

1.06 

5.34 

13.81 

BM (4+4) 

4.28 

17.25 

16.79 


In the process of selecting membrane modes for the basis, the lowest anti-symmetric modes were 
included without regard for the number of periods in the x and y-directions. Since the y- 
dimension is shorter, it is also stiffen This means that membrane modes corresponding to this 
direction will generally have higher natural frequencies than the ones associated with membrane 
motion in the longer plate dimension. Of the first 12 membrane modes used, the maximum 
number of periods in the x-direction was four, while the maximum number in the v-direction was 
two. Consequently, it is expected that modal reduction utilizing 12 membrane modes for the 
plate is less accurate than the best modal reduction for the beam utilizing 6 membrane modes, as 
the latter had six periods along the beam’s span. The conclusion is that the reduced-order 
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analysis for a two-dimensional structure requires an expanded modal basis relative to a one- 
dimensional structure of comparable characteristic length and excitation bandwidth. 


4.2. RANDOM ACOUSTIC RESPONSE 

For the random analysis, a uniformly distributed acoustic load with a bandwidth of 1024 Hz was 
applied. Two loading levels were considered; a low overall sound pressure level of 106 dB 
(5.801 xlO 4 psi RMS) and a high level of 164 dB (0.297 psi RMS). Similar to the beam studies, 
a mass proportional 2% critical damping corresponding to the fundamental mode was prescribed. 
The displacement error corresponding to the random plate response in an RMS sense, as per 
Equation (1 1), is presented in Table 7 for the two loadings levels. 


Table 7: Plate random displacement errors at (%, %) location, %. 



w 

u 

V 

106 dB 

BM (12+12) 

-0.32 

-10.17 

-2.55 

BM (8+8) 

-0.32 

-10.13 

-2.57 

BM (4+4) 

-0.40 

-10.97 

-2.56 

160 dB 

BM (12+12) 

4.00 

-38.68 

-26.43 

BM (8+8) 

14.63 

46.35 

48.02 

BM (4+4) 

20.64 

63.04 

52.00 


For the 106 dB level, which results in a linear response regime with w mm / h = 0.024, the effect of 
modal truncation is negligible; 8, 16, and 24 basis vectors analyses yield essentially the same 
RMS error. As the excitation level is increased to 160 dB, the response becomes highly 
nonlinear with H' max / h = 2.373. At this level, modal truncation starts playing a very important 
role and its effect is significantly magnified when compared to a quasi-static case. The quasi- 
static analysis performed utilizing 8 basis vectors (4+4) for the 0.297 psi load yielded a 
transverse displacement error of about 4%, while the random response analysis for the 160 dB 
load (0.297 psi RMS) required 24 basis vectors (12+12) to obtain comparable results. This trend 
was also observed in the results for the beam. Consistent with observations based on the quasi- 
static analysis, the error in the random transverse displacement is substantially less than the error 
in either component of the random membrane displacement. 
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5. CONCLUSIONS 


Nonlinear modal simulation using an indirect nonlinear stiffness evaluation method has been 
shown to provide accurate predictions of nonlinear quasi-static and random response, when an 
appropriate basis is selected. When used in conjunction with direct numerical simulation for 
validation, the approach constitutes one of the few high-fidelity design options for nonlinear 
random vibration and high-cycle fatigue. 

The following conclusions are limited to the class of problems considered, i.e. those exhibiting 
nonlinear bending-membrane coupling. Of the four modal basis variants considered, the bending 
and membrane modal basis was found to be the only one to accurately predict transverse and 
membrane displacement, and bending and membrane stress at any location on the structure. Its 
main drawback is that identification and selection of membrane modes is labor intensive even for 
the simple beam and plate structures considered. 

The bending-only basis was the simplest amongst the variants in terms of the modal basis 
selection. However, the bending-only basis was incapable of computing membrane 
displacement response, and the accuracy of the membrane stress prediction was highly 
dependent on location. Further, it over-predicted the effects of nonlinearity on the random 
response. 

For the combined bending-companion modal basis, there was no issue with regard to identifying 
the companion modes, and their inclusion did not increase the system size. However, there was 
no engineering rationale for selecting which initial imperfection shape to apply when computing 
the dynamic companion. It was found that results obtained using this basis were comparable to 
those obtained using the bending-only basis, except for the membrane displacement, which the 
latter was unable to predict. The random membrane displacement response did not exhibit 
period doubling, and was significantly reduced in amplitude relative to the physical DoFs 
simulation. 

For the uncoupled bending-companion modal basis, there was also no engineering rationale for 
selecting which initial imperfection shape to apply when computing the dynamic companion. 
When computing the static companion, there was no engineering rationale for selecting which 
displacement or loading to apply. Both static and dynamic UBC bases indicated different 
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coupling between bending and companion modes than the coupling between bending and 
membrane modes for the BM basis. For both static and dynamic UBC bases, a numerical 
instability problem limited the number of companions that could be included in the basis, 
without the use of static condensation in modal coordinates. Finally, for the UBC dynamic basis, 
the membrane displacement and stress response prediction was highly dependent upon location. 

Finally, with regard to two-dimensional structures, it was found that the number of basis vectors 
needed for an accurate response prediction in a highly nonlinear regime had to be considerably 
expanded as compared to one-dimensional structures. 
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